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ABSTRACT 


This thesis presents a study of the atmospheric refractivity conditions 
in the oceanic area around Central America and a description of the possible 
propagation effects on surveillance systems. The basic principles of atmospheric 
refraction are presented along with the techniques used to determine the 
occurrence of atmospheric ducts and to characterize the type of refractive profiles 
occurring at a station. Radiosonde data collected by coastal stations in the area 
were used to establish refractive conditions. 
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I. INTRODUCTION 


A. PURPOSE 

Modern warfare tactics rely heavily on the performance of radar and 
communications systems, both ground based and mobile. Although new 
technology has made significant hardware improvements that have extended 
range detection (power increases, antenna design, digital signal processing, e.g.), 
the atmospheric effects on these systems still play an important role in their 
performance. 

The study of refraction effects has been concentrated in those areas where 
the U.S. military has spent a lot of time, such as the Persian Gulf, the California 
coast and the Mediterranean sea, or areas that present obvious abnormal 
refraction conditions. The purpose of this work is to establish overall refractive 
conditions, determine the availability of data for the region and evaluate the role 
of the near surface refraction effects on the performance of surveillance and ESM 
(Electronic Support Measures) systems for the Central America region. 

B. BACKGROUND 

Many studies have been performed in different areas of the world which 
have resulted in well documented anomalous propagation, especially in the 
Mediterranean, Indian and Northern Pacific oceans. However, an in depth study 
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of refractive conditions does not appear to exist for the Central America's two 
ocean influence region. 

A general meteorological description of the area exists and is given by the 
Forecasters Handbook for Central America and Adjacent Waters (Naval 
Environmental Prediction Facility, Technical Report TR 89-09,1989). However, the 
description has not considered the refraction phenomena of the region and their 
effects on radio-wave propagation in depth. 

Military operations in open ocean and coastal environments depend on 
accurate atmospheric data to estimate systems performance and have to rely on 
the existence of sufficient knowledge of the environment. Unfortunately, data has 
been and continues to be insufficient to meet most planning as well as operational 
needs. The purpose of this thesis is to provide atmospheric information and 
analysis of the Central America region to aid military commands. 

There are computer based models that have been in use for several years 
to assess the effects of the various atmospheric conditions on the performance of 
specific systems. One model is EREPS (Engineers's Refractive Effects Prediction 
System), similar to IREPS (Integrated Refractive Effects Prediction System), which 
is used by the U.S. Navy for operational propagation assessment. Such 
microcomputers models provide a convenient way to investigate the effects of 
particular atmospheric conditions on the performance of systems. 
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II. ATMOSPHERIC REFRACTION 

A. WAVE PROPAGATION 

The propagation of the electromagnetic energy is strongly affected by the 
atmosphere in several ways: refraction, absorption and noise generation 
[Ref. 1] as well as by other phenomena. The latter includes scintillation and 
fading, due to atmospheric inhomogeneities, and turbulence and optical 
interference, due to multi-path propagation. 

There are at least four well differentiated propagation categories: ground 
wave, line of sight, ionospheric and scatter propagation. For low frequency ranges 
(up to perhaps 1500 KHz) ground wave propagation is the main propagation 
mechanism. The ground wave is radiation directly affected by the surface of the 
earth at the interface. In general, it is vertically polarized and very stable over 
time. For higher frequencies, the ionospheric propagation becomes important up 
to a critical frequency near the end of the HF range. Solar induced ionization of 
the upper atmosphere is the dominant factor in the propagation. There are large 
daily and long term variations in the virtual height and strergth of the reflecting 
layers which affect the frequencies that can propagate and the absorption levels 
between any given points. For frequencies above HF the propagation mechanism 
is mainly space wave or line of sight (LOS) propagation. Here, the main loss 
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factor in the propagation is due to the spherical spreading of the wavefront that 
decreases the signal strength as the inverse of the square of the distance. 

The fourth category is scatter propagation that occurs when inhomogeneities 
(in the troposphere) or ions (in the ionosphere) scatter and reradiate the incident 
energy in all directions. A small amount of this energy will eventually reach the 
receiver located beyond the LOS path. This mechanism can lead to reliable, stable 
microwave links as long as enough power is used to overcome the losses. 



This thesis will focus on microwave and, hence, space wave propagation. 
When microwaves propagate through the atmosphere they are affected mainly by 
the lower atmosphere, troposphere, between sea level and approximately 15,000 
meters. This ic the part of the atmosphere that presents the strongest changes in 
pressure, temperature and humidity that determine the refractivity. The most 
severe effect, due to refractive variation, occurs when a layer of air with a 
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negative refractivity gradient bends the rays back toward the earth trapping the 
energy inside a duct. 

B. ATMOSPHERIC REFRACTIVITY 

Deviations from rectilinear wave propagation are determined by the 
variations in the index of refraction of the atmosphere. The main variation of the 
latter is in the vertical direction and is due to the pressure, humidity, and 
temperature changing with height. Refraction occurs because of speed variations 
along the path of the wave that are described by the index of refraction. 

The radio and radar waves propagate through the lower atmosphere at a 
speed, v, lower than the free-space speed as determined by the electromagnetic 
properties of the medium: 

c 

v= _____ 
yju*k 

where: 

c = speed of light 

ji = magnetic permeability 

k = specific inductive capacity 

The ratio of free-space speed to the actual phase propagation speed in the 
medium is defined as index of refraction, n, with a typical value for air of 1.0003 
at sea level. 
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The Equation for atmospheric refractivity is obtained from the Debye's 
theory of the dielectric constant of gases [Ref. 3]: 

£-1 4 n A , 1 p 2 , 1 

e+2 3 M 3 KT V 1 + »cot 


where: 

£ = Dielectric constant = n 2 
A = Avogadro's number 
K = Boltzman's constant 
M = Molecular weight 
P = Electric dipole moment 
t = Relaxation time 
to = Frequency of exiting field 
a = Polarization 
p = Density 
T = Temperature 

At frequencies below approximately 60 GHz., the relative response time of 


the molecules with respect to the frequency makes the product to *x 
with respect to 1. Since e = n 2 and from the equation of state 
[Ref. 2], the expression then becomes: 


very small 



(n -1)=2tc 


A P r 1 M 2 1 

_*_*[cx+_ 

M RT 3 KT 


The pressure P is the sum of the water vapor pressure, e, and the dry air 
pressure, P d . The molecular weight, M, is also the sum of the dry air, and 
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water vapor M w , components. With contribution of these two components 
identified separately, the index of refraction for moist air is: 


/ ^1 ( P ~ e ) ^1 

m-1) =—*--ia.+- 

v ' M, T d M 

a w 


l*(a +C,I) 


where: 

c _2 jcA c _^ 2 

r R z ~3K 

Assigning values to the constants, leads to the refractivity being conveniently 
expressed with the following equation [Ref. 3]: 


N = (1-m)*10 6 


77.6*P^3.73*10 5 *e 

— * J2 


where: 

N = Refractivity 
n = Index of refraction 
P = Pressure, mbar 
T = Absolute temperature, Kelvin 
e = Partial pressure of water vapor, mbar 

The physical composition of the air through the entire troposphere is almost 
constant as far as the oxygen, nitrogen and other gases are concerned. The 
refractivity depends, then, mainly on the temperature and water vapor contained 
in the air. The contribution of the water-vapor term to the refractivity is relatively 
small for cold air because then the saturated vapor pressure is small. But for 
warm air the term becomes important and the water-vapor content has a strong 


influence on the value of N. 






There are no special absorption effects up to 60 GHz., where the oxygen 
molecule shows its first resonance [Ref. 4], well above the frequencies of 
most surveillance radars. 

The water content of the air decreases rapidly with height, as does pressure, 
while the temperature decreases slowly to the tropopause where it stabilizes. 
These changes in the meteorological variables lead to refractivity decreasing with 
height, under normal conditions (Standard Atmosphere), at a rate of 39.37 N units 
per Km. A decrease of N with height means that as height increases, the wave 
propagates faster, closer to the free space speed, causing the rays to bend 
downward following Snell's law as shown in Figure 2. 



To use radar for tactical purposes, the exact knowledge of the path followed 
by the waves and their speed are needed to compute position of the target from 
the timing and angle measured. Since radar systems base their position 
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calculations on U.S. or ICAN "Standard Atmosphere" at any given time, when the 
atmospheric conditions differ, the calculated position of target may be wrong and 
it becomes necessary to include refraction effects in the interpretation of radar 
products. 

The abnormal (anomalous) propagation of electromagnetic waves is called 
subrefraction, super-refraction or trapping, according to the extent and direction 
of the deviation of the rays from standard. 

C. TYPES OF REFRACTIVITY PROFILES 

1. Standard Atmosphere 

The International Civil Aviation Organization (ICAO) has defined a 
standard sounding for the atmosphere at non-tropical latitudes that begins with 
a surface (sea-level) temperature of 15 Celsius. Clearly, it is not applicable to 
tropical latitudes where the sea-level temperature is usually 26 Celsius and, 
therefore, a warmer and, importantly, a more humid Marine Atmospheric 
Boundary Layer (MABL) is expected. Jordan (1958) [Ref. 5] has formulated 
a Mean Tropical Atmospheric Sounding that is useful for calculating a base 
profile of refractivity at tropical latitudes (see Figure 3). 
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Figure 3. Modified Refractivity Profile, Mean Tropical Sounding 
and Low Altitude Detail. 


The low altitude detail in Figure 3 shows the refractive profile in the lower 
atmosphere for both the 1CAN standard atmosphere and that of the mean tropical 
sounding by Jordan [Ref. 5 ]. The refractivity gradient at the low altitude is 
greater for the tropical atmosphere which would lead to greater radio/radar 
horizon. 

The distance, d, to the radio horizon is related to the refractivity gradient 
by: 

d = J2*k*a t *h t 

where, h t , is the transmitter height and k is the ratio of the effective earth's 
radius, a t , to the actual radius of the earth, a : 
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1 


a 


i 


1 t dn 
a + dZ 


The ratio, k, was calculated for the first kilometer in the mean tropical 
sounding as k- 1.722 where for the ICAN standard atmosphere its value is 1.333. 
The radio/radar horizon is proportional to J/2 for a given transmitter height. 
Therefore, on the average, the Jordan tropical sounding has 14% larger radio 
horizon. 

2. Superrefractivity and Trapping 

The atmospheric conditions in which superrefractivity and trapping of 
microwave energy occurs always involve stratification of the atmosphere into 
layers of different index of refraction. Trapping layers and resulting ducts can be 
surface based or elevated. Elevated and surface based trapping layers can be 
detected by radiosonde temperature, pressure and humidity measurements and 
by direct measurement of the refractive profile with an airborne refractometer. 

Superrefractive layers are formed when there is a layer in which the 
temperature increases with height (the temperature decreases with altitude with 
a lapse rate of 9.8 Kelvin per Km in a well mixed atmosphere) and/or the water 
vapor sharply decreases with height. Superrefractive layers become trapping 
layers when the refractivity gradient causes the propagating wave to bend 
downward with a radius of curvature smaller than the radius of curvature of the 
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earth. The energy is reflected at the surface again, (surface based duct), or 
refracted back by the base layer, (elevated duct), trapping it inside the duct as 
depicted in Figure 4. 



Figure 4. Propagation of Radar Energy in a Surface Based Duct. 


Figure 5 illustrates how the different duct types relate to the various 
refractive profiles. 
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Figure 5. Refractive Profiles of the Different Duct Types. 


The evaporation duct is due to an abrupt change in relative humidity at the 
water-air interface, from 100 % to approximately 85 % in the first few meters. The 
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evaporation duct is almost always present over water. Its thickness varies as a 
result of the daily temperature changes and intensity of the surface wind. The 
evaporation duct height can be computed from bulk meteorological measurements 
(sea surface temperature, air temperature and wind) by using semi-empirical 
models like the model described by Paulus [Ref. 6]. For a given duct 
thickness, d, the minimum frequency, F^, that can propagate within the duct 
is given by [Ref. 4]: 



l dz) 


The minimum frequency propagated in the surface based duct ?s a function 
of the duct thickness can be seen in Figure 6. 
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D. METEOROLOGICAL PHENOMENA LEADING TO TRAPPING 

Several meteorological conditions modify the distribution of temperature and 
water vapor in the lower atmosphere. Variations in temperature, pressure, or 
both, determine as described before, the atmospheric refractive gradient. Air flow 
associated with the motion of air masses, radiation, subsidence and land-sea 
winds, define the vertical temperature and moisture profile of the atmosphere. 

a. Fadiation 

Radiative cooling of the surface during clear nights is a very 
common cause of temperature inversion leading to superrefractive conditions over 
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land. At sea, the ocean acts as a mass’ 1 ■» heat reservoir so its surface temperature 
remains nearly constant and the effect is diminished. 

b. Subsidence 

Subsidence, the slow, large scale sinking of the air, brings dry air 
from aloft and raises the temperature of the air by adiabatic compression. This 
creates inversion layers that are common on large tropical and subtropical ocean 
areas. Subsidence occurs within subtropical high pressure systems. 

c. Land-sea winds 

The difference in temperature and humidities between the land, 
heated by the sun, and the adjacent sea surface usually causes the warmer, dry 
air over land to rise and subside over the adjacent cooler and more moist marine 
air. The colder moist air from the sea will replace the raising air as depicted in 
Figure 7. The warm, dry air will form the upper boundary of an inversion layer 
in which ducting can occur. 




Sun 

Radiation 



Warm dry air 


* * 


Cold moist air 



Figure 7. Sea-Wind Inversion Layer. 
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III. CENTRAL AMERICA CLIMATOLOGY 


A. CLIMATOLOGY 

The meteorology of Central America is significantly influenced by the 
Intertropical Convergence Zone (ITCZ) that lies on the pacific side most of the 
year [Ref. 7]. The large scale air flow along the ITCZ is the product of the 
northeasterly trade winds meeting the southeasterly trade winds. Since they are 
air masses of different origin, their temperatures and water vapor content are 
quite different making the ITCZ a zone of very high cloudiness and precipitation 
that can be seen all year around in the satellite imagery. Rainfall occurs primarily 
between June and September and by December with drier periods in between. 

The general air flow direction along the ITCZ is west-to-east during most of 
the year. It is deviated by the Andean ranges at the western coast of South 
America, which have heights up to 5000 meters. The Caribbean side is mostly 
affected by the flow of trade winds that flow in harmony with the genera! 
Caribbean regime. Here the dominant phenomena are the tropical disturbances. 
Sea surface temperatures are high in both Pacific and Caribbean areas, with small 
diurnal and annual changes [Ref. 8], which in turn cause very constant 
surface air temperatures. 
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B. SYNOPTIC DESCRIPTION 


The Pacific and Caribbean sides of Central America are subject to different 
regimes. The Pacific presents a general surface air flow following the ITCZ 
parallel to the equator that ends at the western coast of South America where is 
deviated by the high Andes mountains. The Caribbean is farther north and is less 
affected by the ITCZ. 

Local phenomena, such as like sea-breeze and land-sea winds are known to 
cause strong inversion layers in coastal areas. The breezes are frequent and strong 
enough to dominate the general circulation. The complete synoptic description of 
the area can be found in the Forecasters Handbook for Central America and 
Adjacent Waters [Ref. 7]. 
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IV. DATA ANALYSIS AND RESULTS 

A. RADIOSONDE DATA 

1. Geographical Distribution of Data 

For the purpose of this study, we have considered data taken during 
1988 and 1989 by the stations shown in Figure 8. These are listed in Table 1. 



Figure 8. Geographical Location of the Sounding Stations. 
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Table 1 


Stations Used in the Analysis 


Station 


Lat. 

Long. 

Elev. 

78384 

Roberts Field, Gran Cayman Island 

18-19N 

81-20W 

3 m 

78397 

Kingston, Jamaica 

18-04N 

76-50W 

1 m 

78583 

Belize International Airport, Belize 

17-36N 

88-18W 

5 m 

78806 

Howard A.F.B., Panama 

8-58N 

79-35W 

66 m 

80001 

San Andres Island, Colombia 

12-34N 

81-41W 

2m 


There is only one station located on the Pacific side of Central America 
(station 78806 in Howard A.F.B., Panama) in the area of this study. Two of the 
five stations located in the Caribbean are in the Greater Antilles (stations 78384, 
78397 in Jamaica and Gran Cayman), one is in an open ocean island (station 
80001, San Andres) and one is a coastal station in Belize (station 78583). 

2. Distribution of Data Over Time 

Stations are supposed to launch radiosondes twice a day (stations 78806, 
78583, 78397, 78384) and daily in the case of station 80001 [Ref. 9]. Due 
to different circumstances they generally do not. The resultant distribution of 
soundings is presented according to stations in Table 2. 
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B. ANALYSIS OF REFRACTIVE CONDITIONS 


1. Procedure of Analysis 

The Fleet Numerical Oceanography Center, Monterey, CA, provided 
significant and mandatory level radiosonde data on 9-track tape for 1988 and 
1989. The tape files were read and transferred to a PC-type microcomputer for 
processing. Pressure, temperature and dew-point depression from significant 
levels were used to calculate refractivity and height for display as a refractivity 
profile and to compute the refractivity gradient to identify anomalous conditions. 

Since ducts occur at various altitudes, with a broad strength range, it is 
difficult to appropriately describe with statistical measures the refractive profile. 
Therefore, each refractive profile in which ducting was detected was displayed 
and visually evaluated with the purpose of staling a typical ducting profile for 
each station. EREPS was used to graphically show the effects of these profiles on 
a typical radar system performance. 

The ray tracing subprogram of EREPS (RAYS) was used to visualize effects 
of the typical refractive anomalies detected in the soundings. The propagation 
model (PROP) [Ref. 10] was used to assess the performance of typical 
systems under the ducting conditions detected. Finally, the evaporation duct 
climatology from EREPS is used to visualize the expected effect of the evaporation 
duct in surveillance radar and ESM systems. 
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2. Formulas and Algorithms 

Formulae used in the programs to identify and plot ducts are presented 
here and the general procedure followed to determine the occurrence of ducts in 
the soundings. 

Modified Refractivity is calculated from pressure, temperature, partial 
pressure of water vapor and height with the following equation: 

77.6 *P 3.73*10* *e h 

M -- + -+- 

Ta To 1 6.371 

where: 

P = Pressure in millibars 
Ta = Temperature in Kelvin 
e = Partial pressure of water vapor in millibars 
h = Height in meters 

Partial water vapor pressure was calculated from the saturated vapor 

pressure at a the dew point temperature. The saturated water vapor pressure for 

a given temperature is calculated using Lowe's polynomial [Ref. 11]: 

Es = /i 6 *T 6 +/ij*r 5 +/i 4 *7 4 +A3*r 3 +A 2 *r 2 +A*r+i4 0 


where: 

A 0 = 6.107799961000000 
A = 0.443651852100000 
A 2 = 0.014289458050000 
A 2 = 0.000265064847100 
j4 4 = 0.000003031240396 
A s = 2.0340809481 *lO* 

A 6 = 6.136820928999999 * 10 n 
T = Temperature in Celsius 
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The Height is calculated from the initial launch altitude and pressure for the 
first level, and from the preceding level for the successive heights using the 
relationship for geopotential height [Ref. 12], as expressed in the 
following equation: 

h = V14-63 

“i 


where: 

h 0 = Initial Height, meters 

Vt 0 = Virtual Temperature at Initial Height, Kelvin 
P 0 = Pressure at Initial Height, mbar 
Vt x = Virtual Temperature at Actual Height, Kelvin 
Pj = Pressure at Actual Height, mbar 

Virtual temperature is computed from the temperature, pressure and water 

vapor pressure with the equation: 


Vt _ Ta , (0.3794017*7fl*e) 
iP-e) 


where: 

Ta = Absolute temperature (Kelvin) 
P = Pressure 
e = Vapor pressure. 


The trapping layers are identified by computing the refractive gradient 
dM/dZ and plotting the refractive profile, modified refractivity, M, versus height, 


25 









r" 


for the data sets that have negative gradient to visually evaluate the resulting 
trapping layers. 

3. Results by Stations 

a. San Andres Island, Station 80001 

A total of442 soundings where processed for this station. Trapping 
layers were classified as being surface based or not, and if not, as being above or 
below 700 mb. Trapping layers were found in 156 soundings, and are distributed 
as follows according to occurrence being at the surface or, if not, as being above 
or below 700 mb : 

l 


Year 

1988 

1989 

Total 

Soundings 

161 

281 

442 

Surface based trapping layers 

15 

16 

31 

Elevated trapping layers (total) 

41 

84 

125 

Elevated trapping layers above 700 mb. 

10 

20 

30 

Total 

56 

100 

156 


The refractive profiles that showed trapping layers were visually evaluated 
and it was found that trapping layers were stratified much like the profile in 
Figure 11. Many profiles showed more than one trapping layer and therefore have 
multiple ducts. Only the lower trapping layer per sounding was used for the 
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above statistics. The trapping layer height was about 2000 meters and below. The 
thickness and intensity of trapping layers were less than or equal to the profiles 
shown. The local launch time for this station is 06:30 (12:00 UT). Therefore the 
observations arc representative of a stable, non mixed MABL in which 
stratification is enhanced. 

Soundings revealing trapping layers had many data errors. Nineteen of the 
156 layers identified by the program were due to errors in the data. These were 
identified visually by checking the suspect part of the data for robustness. Surface 
based trapping layers were found to occur in 7 % of the soundings. Because 
surface trapping layers could have been caused by erroneous launch level values, 
they were checked closely. The total duct occurrence found is 31 %, after 
excluding the errors. For this station, Ortenburguer and Lawson [Ref. 13] 
has found an occurrence of surface based ducts of 20 % and an occurrence of 


elevated ducts of 35%. 









80001 — Radiosonde Data 1988 



Figure 9. Radiosonde Data, 1988, San Andres, Colombia, Station 
80001. 


80001 — Radiosonde Data 1989 



1 2 3 4 5 6 7 8 9 10 11 12 


Month 

Figure 10. Radiosonde Data, 1989, San Andres, Colombia, Station 
80001. 
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Figure 11. Refractivity Profile, San Andres Island, 
Station 80001, 25 Jan 1988. 
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Figure 13. Refractivity Profile, San Andres Island, 
Station 80001, 14 Mar. 1989. 


b. Howard A.F.B., Panama, Station 78806 

For this station 766 soundings where processed, trapping layers 
were found in 266 soundings distributed as follows: 


Year 

1988 

1989 

Total 

Soundings 

417 

359 

776 

Surface based trapping layers 

45 

45 

90 

Elevated trapping layers (total) 

107 

69 

176 

Elevated trapping layers above 700 mb. 

25 

16 

41 

Total 

152 

114 

266 
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The soundings that showed trapping layers were visually evaluated. The 
occurrence of trapping layers was observed centered at around 2500 meters and 
the associated ducts tend to be strong and thick as seen in the profile of Figure 16. 
The strength of the duct refers to the M difference across the duct while thickness 
refers to the vertical separation of the duct top and base. Many refractive profiles 
contain more than one trapping layer and associated duct. Local launch times are 
06:41 and 18:41 (12:00 UT and 00:00 UT). Data from these times could have 
greater ducting occurrence than mid-day time data due to lesser convection 
generally occurring before sunrise and after sunset. 

Visual inspection of the first 5000 meters in the profiles that showed 
trapping layers, revealed fifteen data sets in error of the total 266 soundings with 
trapping layers found by the program. This reduces the duct occurrence found 
from 34.7 % to 32.3 %. The surface based inversion layer occurred in 11.6 % of 
the soundings. As was described above, a surface based inversion layer could 
arise from erroneous surface values. This possibility was checked very closely for 
this case because the percent occurrence of surface based trapping layers was 
above 10%. For this station, Ortenburguer and Lawson [Ref. 13 ] has found an 
occurrence of surface based ducts of 5 % and an occurrence of elevated ducts of 
29 %. 
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78806 - Radiosonde Data 1988 



Figure 14. Radiosonde Data, 1988, Howard AFB., Panama, Station 
78806. 


78806 — Radiosonde Data 1989 



Month 


Figure 15. Radiosonde Data, 1989, Howard AFB, Panama, Station 
78806. 

























Figure 16. Refractivity Profile, Howard AFB, 


Panama, Station 80001, 25 Mar. 1989. 



Figure 17. Refractivity Profile, Howard AFB, 


Panama, Station 78806, 13 Apr. 1988. 
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c. Belize Airport , Belize; Station 78583 

A total of 615 soundings where processed for this station. Trapping 
layers and associated ducts were found in 188 soundings and are distributed as 
follows: 


Year 

1988 

1989 

Total 

Soundings 

397 

218 

615 

Surface based trapping layers 

10 

7 

17 

Elevated trapping layers (total) 

108 

63 

171 

Elevated trapping layers 700 mb. 

36 

13 

49 

Total 

118 

70 

188 


The elevated trapping layers occurred consistently at an altitude of about 
1700 to 2000 meters. There where few multiple duct occurrences and the few 
surface based ducts where not very strong. Figure 20 and Figure 21 show typical 
M profiles. The inversion layer is stronger and thicker in the 12:00 UT sounding. 

Two errors were detected in the data of the soundings found with ducts. 
The percentage of trapping layer occurrence was 30.2 %. Surface based inversion 
layers occurred in 2.8 % of the soundings. This station is not included in the work 
by Ortenburguer and Lawson [Ref. 13]. 
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Month 


Figure 18. Radiosonde Data, 1988, Belize Airport, Station 78583. 


78583 — Radiosonde Data 1989 



Month 


Figure 19. Radiosonde Data, 1989, Belize Airport, Station 78583. 




















Figure 20. Refractivity Profile, Belize, Station 
78583, 14 Mar. 1989. 
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Figure 21. Refractivity Profile, Belize, Station 
78583,14 Mar. 1989. 
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d. Kingston, Jamaica, Station 78397 


A total of 1182 soundings were 

processed for this station. Trapping 

layers and associated ducts were found in 248 soundings and 

follows: 

are distributed as 

Year 

1988 

1989 

Total 

Soundings 

659 

523 

1182 

Surface based trapping layers 

24 

16 

40 

Elevated trapping layers (total) 

123 

85 

208 

Elevated trapping layers above 700 mb. 

25 

27 

52 

Total 

147 

101 

248 


Visual evaluation of the 147 refractive profiles that showed trapping layers 
revealed nineteen erroneous data sets. The trapping layer is located generally at 
2000 to 2500 meters. Most ducts are neither strong nor thick. The profiles shown 
in Figure 24 and Figure 25 are representative of the strongest trapping layers seen. 
The frequency of occurrence of trapping layers was 19.4 %. Surface based 
trapping layers occurred in 3.4 % of the soundings. For this station, Ortenburguer 
and Lawson [Ref. 13] has found an occurrence of surface based ducts of 20 % 
and an occurrence of elevated ducts of 33 %. 
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78397 — Radiosonde Data 1988 
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Figure 22. Radiosonde Data, 1988, Kingston Jamaica, Station 78397. 


7839 7 - Radiosonde Data 1989 
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Figure 23. Radiosonde Data, 1989, Kingston, Jamaica, Station 78397. 
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Figure 24. Refractivity Profile, Kingston, Jamaica, 
Station 78397, 5 Jan. 1988. 
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e. Roberts Field, Gran Cayman, Station 78384 


For this station, 1148 soundings where processed and 381 
soundings were found to contain trapping layers distributed as follows: 


Year 

1988 

1989 

Total 

Soundings 

618 

530 

1148 

Surface based trapping layers 

19 

18 

37 

Elevated trapping layers (total) 

194 

150 

344 

Elevated trapping layers above 700 mb. 

47 

25 

72 

Total 

213 

168 

381 


The visual evaluation of the refractive profiles that contained trapping layers 
revealed that the elevated trapping layers where thick and strong and were 
grouped mainly near 1500 or near 3000 meters. 

The occurrence of trapping layers was found to be 33.2 % and the 
occurrence of surface based trapping layers was 3.2 %. 



78384 — Radiosonde Data 1989 
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Figure 27. Radiosonde Data, 1989, Gran Cayman Island, Station 
78384. 
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Figure 28. Refractivity Profile, Gran Cayman, Sta 
tion 78384,14 Mar. 1989. 
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Figure 29. Refractivity Profile, Gran Cayman, 
Station 78384,13 Die. 1988. 








V. PROPAGATION EFFECTS ON SYSTEMS 


A. EFFECTS ON RADARS 

The radar equation that conveniently describes the performance of a radar 
in terms of transmitted power, antenna gain and effective aperture, radar cross 
section of target and receiver sensitivity falls short when used to describe actual 
performance of systems. 

"In practice, however the simple radar equation does not predict the 
range performance of actual radar equipments to a satisfactory degree of 
accuracy. The predicted values of radar range are usually optimistic. In 
some cases the actual range might be only half that predicted." [Ref. 4] 

The radar equation describes the free space range and assumes "...three- 
dimensionality and isotropy of space" [Ref. 1], However, the atmosphere is in 
no way isotropic and the three-dimensionality is violated as soon as ducting or 
other anisotropic phenomena are encountered. Super refraction and ducting can 
direct energy in a particular direction that otherwise would spread 
isotropically producing extended coverage, sometimes to ranges well beyond 
free-space ranges. 

Extended range is not the only effect expected. Due to the reciprocal path, 
there is a significant increase in the clutter and noise received that can 
seriously degrade the overall performance of a surveillance system. Double 
return echoes are possible and a target at greater than the maximum 
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appear as an in-range target. The distance at which the signals can be intercepted 
will be augmented and the scenario becomes much more complex. 

The surface-to-surface propagation over the ocean is affected in the area due 
to the moist MABL. The radio horizon (Page 10) is slightly extended on the mean 
over what would be expected over subtropical waters and somewhat larger than 
the world average. 

1. Rain Effects 

Rain is the most common and an important phenomena in tropical areas 
as nearly all tropical disturbances carry rain. Absorption and scattering of the 
radar energy by rain is a very strong loss mechanism. The difficulties involved 
in accurately measuring the characteristics and the rate of rainfall prevent the 
models from being effectively used to forecast the effect of rain on the radar 
coverage. MTI (Moving Target Indicator) or doppler radars show less rain clutter 
effects; but, since the energy is scattered in the rain volume, the effective power 
that reaches the target and the receiver is strongly reduced and, therefore, so is 
the detection range. The rain scatter is frequency dependent [Ref. 14] so 
low frequency radars can be used if see-through capability is needed. 

The transient nature of rain, even in high precipitation areas like the Pacific 
ocean off the western coast of Colombia, suggest the selective use of sensors as 
the best way to obtain radar coverage during rain. 
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The transient nature of rain, even in high precipitation areas like the 
Pacific ocean off the western coast of Colombia, suggest the selective use of 
sensors as the best way to obtain radar coverage during rain. 

2. Duct Effects 

The range of radar coverage can be extended by surface based ducts 
associated with elevated trapping layers or by evaporation ducts. The 
evaporation duct is a subset of ducts where the duct base is at the surface. In 
the a 1 sence of a surface based duct, the evaporation duct is tb° dominant 
propagation mechanism for surface-to-surface detection over water, provided 
the radar antenna and the target are inside the duct and a frequency 
sufficiently high is used (above the cutoff frequency). 

a. Coverage Diagram with Surface Based Duct 

Due to the large number of different ducting profiles observed 
at the stations analyzed, the effects a surface based duct on a particular 
system ducts will be illustrated using a sounding at station 80001 on March 14 
1989. The refractive profile obtained from the sounding is depicted in 
Figure 13. A strong trapping layer is clearly visible from the surface up to 100 
meters, causing a surface based duct. At around 3000 meters there is another 
thin duct. As stated earlier f^r this station, many soundings in which a duct is 
found revealed multiple ducts. 
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inside the duct. For comparison purposes, a standard refraction raytracing is 
shown in Figure 31. 
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Figure 31. EREPS Ray Tracing of a Standard Refractive 
Atmosphere. 


b. Propagation Loss 

The ray tracing alone does not give enough information on the 
detection capability of a radar under specific refractive conditions. The system 
itself needs to be taken into account. Receiver characteristics as well u3 antenna 
gain, transmitter power, target size, thresholds and other factors need to be 
considered to assess detection capability. Table 3 lists the parameters necessary 
to estimate loss. The end application of the radar will ultimately dictate the 
probability of detection and false alarm rate to be included in the model. 

Here, for the purpose of visualizing the effects of ducting in a given system, 
a "generic" surveillance radar is used with a target of small radar cross section 


47 





(RCS) and a Swerling model 1 [Ref. 4] to account for the slow fluctuation of RCS 
typical of a low tonnage vessel at sea. 


Table 3 

Parameters Used to Estimate Propagation Loss 


Frequency 

3 GHz. 

Polarization 

Horizontal 

Radar height 

10 meters 

Target height 

15 meters 

Antenna type 

Cosec-Sq. 

Vertical BW. 

10 degrees 

Elevation Angle 

0 degrees 

Peak power 

500 Kw. 

Pulse width 

1.0 uSec. 

Antenna gain 

32 dBi. 

System losses 

8.4 dB 

Receiver noise figure 

5dB 

Pulses integrated 

10 

RCS target 

10 sqm. 

Probability of detection 

.8 

Probability of false alarm 

io- 8 

Swerling case 

1 


The effects of the duct, shown in Figure 32, are tactically important. There 
is a significant skip zone in which there is no detection. The detection range is 
extended beyond the free-space range after the first reflection at 30 nm up to 55 
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nm. This an additional coverage of 6000 squared nautical miles to the 346 nm. 2 
of the unextended coverage. If a force is aware of the presence of the duct, it can 
be used on purpose to get extended surface surveillance. Otherwise, the force can 
give away valuable intelligence by ignoring the existence of the duct. 



Figure 32. Propagation Loss Plot for the Duct found on San 
Andres Island, Station 80001 on 3 Jun. 1989. 


The Evaporation duct exists nearly all the time. The summary data provided 


by EREPS calculated from 15 years of observations shows that an evaporation 


duct height of more than six meters 94.5 % of the time (Figure 33) and 16 or more 


meters high 53.1 % of the time. These are high values compared with the world 
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average (Figure 34). The evaporation duct is an important propagation 
mechanism, particularly for higher frequencies. Tactical effects of the 
evaporation ducts are often neglected and underestimated. Since the duct is 
present most of the time and is both frequency and height dependent, small 
vessels with their radar located inside the duct, operating in the X band, may 
have earlier detection of conventional platforms with higher antennas and 
lower frequency radars. 





Figure 33. Evaporation Duct Summary, San Andres Island, 


Station 80001. 



Figure 34. Evaporation Duct Height, World Average. 
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The range of the same radar in the presence of an evaporation duct of 16 
meters is extended from 10 nmi. to 15 nmi., an increase in coverage from 314 
square nmi to 706 square nmi. The increase in area coverage (225 %) and the 
relative permanence of the duct (expected 53% of the time at station 80001) 
make it an important propagation mechanism. The effects are shown in 
Figure 35. 



Figure 35. Extended Coverage Due to Evaporation Duct. 


B. EFFECTS ON ESM 


Anomalous propagation effects on ESM can be dramatic. Anomalous 


propagation can result in signal strength that are orders of magnitude bigger 


or smaller than the free-space signal, rendering relative signal strengths useless 
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as a rough distance measure; and the attenuation effects are not readily 
observed in the displays due to the passive nature of the system. 

Today's expected operating environment is so dense that digital 
processing techniques are always used in conjunction with ESM and radar 
warning receivers (RWR) to sort and extract the various parameters of a signal. 
Sensitive ESM receivers suitable for electronic intelligence gathering (ELINT) 
will eventually operate in such a high densities that the processing cannot cope 
with the task. The probability of intercept for a given signal will be greatly 
reduced. 

The saturation of the processing can result from any mechanism that 
effectively brings more signals into the receiver. Altitude, in an airborne 
receiver, and ducts, for both surface and air borne systems, effectively extend 
the radius of coverage for a given sensitivity. Assuming a given number of 
emitters per unit area, the number of emitters received will follow the square 
of the range covered, resulting in a substantial increase of the density of the 
environment. What could be seen as an advantage of getting extended 
coverage can very well result in a lower probability of intercept for all signals. 






VI. CONCLUSIONS 


The atmospheric effects on VHF/UHF/Microwave propagation in the 
tropical waters surrounding Central America are significant. There is a high 
probability of ducting and extended radar horizon is expected most of the 
time. 

Statistical and historical data are based on the soundings taken by the 
WMO station network. These soundings are used for the synoptic scale 
weather analysis and are not intended for local refractivity study. Data 
collection is seriously flawed when the times, instruments and collection 
procedures are not precise and consistent enough to draw conclusions on the 
more dynamic refractivity behavior. There is no substitute for real time, 
tactically planned soundings and data collection. 

The data available for the area are therefore just a sample taken at the 
few stations shown and at particular times in which ducting is more probable. 
The analysis of the results of refractivity studies using these data should be 
done using extreme caution. The overall tropical environment must be 
evaluated for consistency. Routine software checks are not enough when 
dealing with atmospheric data. Knowledge of the physics involved in the 
atmospheric behavior is necessary for military personnel due to the operational 
consequences of incorrectly assessing refractive conditions. 
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Further studies are needed to validate the interpretation of the refractive 
data available for the area against measurements or case observations. There 
are few seasonal variations in the temperature and wind in the tropical area 
^ and therefore correlation with satellite imagery may be useful in determining 

the extent of anomalous refractive conditions. 

There is the need for increasing the awareness of the military on the 
effects of propagation on the systems and tactics to be used, as well as the 
particular climatology of the operating theater in Central America. The study 
of the environment and its effects on radio propagation should be included in 
the undergraduate education and intensified in the graduate education to 
provide a wide base of knowledge and understanding of systems capabilities. 

The only radiosonde data available for the Pacific Ocean comes from the 
4 Howard Air Force Base, Panama Canal Zone, which, as seen in the map in 

Figure 8, cannot be very representative of the conditions in Central America's 
Pacific Ocean area. Not enough data is available to even characterize the 
refractive phenomena that may occur there. Alternative sources such as 
satellite imagery and remote sensing of the atmosphere could be used to 
obtain the data that current environmental models use. 
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APPENDIX A. 


f 


The programs to manipulate the radiosonde data were written in Quick- 
Basic^) version 4.5 and are too big to list, the routines used to obtain 
modified refractivity and height and the general subprogram used to identify 
and plot ducts, are presented here. 


' P Pressure in millibar 

' Ta Temperature in Kelvin 

' e Partial pressure of water vapor in millibar 

' H Height in meters 

/ 

FUNCTION RefractM (p, Ta, e, H) 

M = 77.6 * p / Ta + 3.73 * 10 A 5 * e / (Ta * Ta) 
+ H / 6.371 

RefractM = M 

END FUNCTION 


' Function to find the saturated vapor pressure (millibar) 
' using Lowe's polynomial 

' Temperatures in Celsius 
/ 

FUNCTION Satvp (Temp) 

DEFDBL A 

AO = 6.107799961# 

A1 = .4436518521# 

A2 = .01428945805# 

A3 = .0002650648471# 

A4 = .000003031240396# 
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A5 = 2.0340809481D-08 
A6 = 6.136820928999999D-11 

Svp = (((Temp * A6) + A5) * Temp + A4) 

Svp = ((Svp * Temp + A3) * Temp + A2) * Temp + A1 
Svp = Svp * Temp + AO 
Satvp = Svp 

END FUNCTION 


' HO Initial Height, meters 

' VtO Virtual Temperature at Initial Height, Kelvin 
' PO Pressure at Initial Height , mbar 
' Vtl Virtual Temperature at Actual Height, Kelvin 

' PI Pressure at Actual Height , mbar 
/ 

FUNCTION Height (HO, VtO, Vtl, PO, PI) 

Height = HO + 14.63 * (VtO + Vtl) * LOG(PO / PI) 
END FUNCTION 


' Ta Absolute temperature (Kelvin) 

' PI Pressure 

' Ee Vapor pressure environ. 

/ 

FUNCTION Tv (Ta, PI, Ee) 

Tv = Ta + (.3794017# * Ta * Ee) / (PI - Ee) 
END FUNCTION 


/ 

SUB Process (Alt$, Layer$, Intensity, DuctHeight, NrData%) 
Layer $ = " " 
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Tmp$ = "" 

Intensity = 0 
LayerHeight = 0 
PO = ppp(O) 

HO = VAL(Alt$) 

H(0) = HO 

FOR i = 0 TO NrData% - 1 
Td = TTT(i) - TTD(i) 

Ee = Satvp(Td) 

Ta = AbsoluteTemp(TTT(i)} 

PI = ppp(i) 

Vtl = Tv(Ta, PI, Ee) 

IF i = 1 THEN 
VtO = Vtl 
END IF 

H(i) = Height(HO, VtO, Vtl, PO, PI) 

M(i) = RefractM(PI, Ta, Ee, H(i)) 

HO = H (i) 

PO = PI 
VtO = Vtl 

IF H(i + 1) - H(i) = 0 THEN 

Duct$ = "X" 
dMdZ(i) = 0 
GOTO OnZero 
END IF 

dMdZ(i) = (M(i + 1) - M(i)) / (H(i + 1) - H(i)) 

OnZero: 

IF dMdZ(i) < 0 THEN 

SELECT CASE ppp(i) 

CASE IS > 700 

Layer$ = "D" 

LayerHeight = H(i) 

IF i = 0 THEN 

Layer$ = "S" 

END IF 

CASE IS <7 00 

IF Layer$ = " " THEN 
Layer $ = "d" 

LayerHeight = H(i) 

END IF 
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END SELECT 
END IF 

IF dMdZ(i) < Intensity THEN 
Intensity = dMdZ(i) 

END IF 
NEXT i 

'Plot the Layers found below 700 mbars. 

IF Layer$ = "D" OR Layer$ ■ "S” THEN 
Plot M(), H(), NrData%, Layer$ 
END IF 
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